Three-hadron angular correlations from pQCD at RHIC and LHC 
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We study three-hadron azimuthal angular correlations in high energy proton-proton and central 
nucleus-nucleus collisions at RHIC and LHC at mid-rapidity. We use the LO parton matrix elements 
for 2 — j> 3 processes and include the effect of parton energy loss in the Quark-Gluon Plasma using 
the modified fragmentation function approach. For the case when the produced hadrons have either 
same or not too different momenta, we observe two away side peaks at 2tt/3 and 47r/3. We consider 
the dependence of the angular correlations on energy loss parameters that have been used in studies 
of single inclusive hadron production at RHIC. Our results on the angular dependence of the cross 
section agree well with preliminary data by the PHENIX collaboration. We comment on the possible 
contribution of 2 — >■ 3 processes to di-hadron angular correlations and how a comparison of the two 
processes may help characterize the plasma further. 
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PACS numbers: 25.75.-q, 25.75.Gz, 12.38.Bx 



I. INTRODUCTION 

Jet quenching, the energy loss of a parton moving in a 
Quark-Gluon Plasma (QGP), is one of the key paradigms 
emerging from high energy heavy ion collisions at the 
Rclativistic Heavy Ion Collider (RHIG) and the Large 
Hadron GoUider (LHG). This energy loss is believed to 
cause the suppression of single hadron transverse momen- 
tum spectra in Au -f- Au collisions, as compared to p -f- p 
collisions, as well as the disappearance of the away-side 
peak in two-particle azimuthal correlations [l| . Neverthe- 
less, the fact that many different energy loss models can 
fit the observed suppression of single and double inclusive 
hadrons, produced in high energy heavy ion collisions, 
makes it imperative to consider more exclusive observ- 
ables in order to shed light on the dynamics of energy 
loss. One can ask for example whether the double hump 
structure on the away-side in azimuthal two-particle cor- 
relations in Au -I- Au collisions [2| can be studied from 
looking at more exclusive channels such as three-hadron 
correlations. These studies have shown to provide a pow- 
erful tool to distinguish between different proposed sce- 
narios 0,0]. 

This question has been addressed in Ref. [5[ where we 
showed that some properties of the two-particle correla- 
tions can be understood from three-hadron production 
studies. The connection is made by considering the dif- 
ferent path lengths traveled by the two away side partons 
in 2 — > 3 versus the away side parton in 2 — >■ 2 parton 
processes. In the former, one of the two final state par- 
tons in the away side travels, on the average, a larger 
path length than the other and also larger than the away 
side parton in the latter processes, leading to an increase 
in the probability of absorption of this parton due to en- 



ergy loss. Since the parton that travels the smaller path 
length in 2 — > 3 processes in the away side has on the 
average, also a smaller path length to go through than 
the away side parton in 2 ^ 2 processes, the parton in 
the former has a smaller chance of loosing energy than 
its counterpart in the latter processes. The combination 
of these two effects effectively amplifies the production 
of structures that upon accumulation of signal look like 
a double hump or a broader peak in the away side, de- 
pending on the momentum difference between leading 
and associate particles. 

To look closer at this issue, in this work we consider 
the three-hadron production cross sections in p -I- p and 
A -I- A collisions using leading order (LO) 2 — )■ 3 par- 
tonic processes. To include the effect of QGP and parton 
energy loss on the final particle spectra, we use modified 
fragmentation functions. We then investigate the effect 
of energy loss on the angular correlations of the three 
produced hadrons by varying the parameters associated 
with the energy loss suffered by each parton. We com- 
pare our results to the preliminary data on three-hadron 
correlation function measured at RHIG and find that the 
angular dependence of the correlation function is repro- 
duced correctly by our formalism. 

The work is organized as follows: in Sec. II we write 
the three-hadron production cross section in p -I- p and 
also introduce the modified fragmentation function which 
takes parton energy loss into account for A -I- A collisions. 
We show the time-evolved profile of the medium and the 
average number of scatterings as functions of the initial 
points and trajectories on the transverse plane for central 
collisions. In Sec. HI we investigate the dependence of 
the three-hadron production cross section at RHIG and 
LHG on energy loss parameters. We then compare our 



results to preliminary data on the angular correlations 
of the three produced hadrons and show that there is 
a good agreement in the kinematic range accessible by 
our formalism. We conclude in Sec. IV by pointing out 
to the possible contributions of 2 — > 3 processes to di- 
hadron correlations and how they may probe different 
transverse plane evolution of the hard scatterings within 
the plasma. 



II. THREE-HADRON PRODUCTION 



where y/s is the total center of mass energy available for 
the collision, xi is the momentum fraction of the incom- 
ing parton in the projectile. The momentum fraction of 
the incoming parton in the target is fixed also by mo- 
mentum conservation to be xi. The angular dependent 
part of the phase space factor is 
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The differential cross section for three-hadron produc- 
tion, with momenta hi, /12, /13, in p -I- p collisions in 
mid-rapidity is given by 
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where /j / and /,- /p are the distribution functions of par- 
tons i, j within the colliding protons. We use the CTEQ6 
parametrization [6|. Aiij^p^p^p^ is the corresponding 
leading order matrix element describing the process at 
the parton level Q- 

The sum is over all the possible colliding parton 
species, over all of the three possible partons in the final 
state and over all the permutations for a given parton in 
the final state to become a given hadron. The magni- 
tudes of the final state parton momenta are determined 
from momentum conservation at the parton level as func- 
tions of the angles that p2 and ps {62, 03, respectively) 
make with pi 
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and Z?p^ /jj^^ (/i,„/p„), is the fragmentation function of 
parton P„ to become hadron H„i, which is a function 
of the momentum fraction Znm = hm/pn {n,m = 1,2, 3). 
We take the fragmentation functions for charged hadrons 
as given by the KKP parametrization [8| . In each permu- 
tation, any given parton can become the leading hadron 
which, wc take as Hi, and the other two the away-side 
hadrons. We work in the limit of collinear fragmentation 
thus, the angles that define the direction of the away- 
side hadrons, 0f {i = 2,3), are linearly related to the the 
parton angles 0j {j = 2,3). 

To consider the process within a central heavy-ion col- 
lision and thus account for the effects of energy loss, we 
resort to the model put forward in Ref. Q- The model 
considers an initial gluon density obtained from the over- 
lap of two colliding nuclei, each with a Woods-Saxon den- 
sity profile. The gluon density of the medium is diluted 
only due to longitudinal expansion of the plasma since 
transverse expansion is neglected. 

The gluon density pg is related to the nuclear geometry 
of the produced medium as 



Pg{T,b,r,n) 



T 2A 

TA{\r + hT\)+TA{\b-f-nT\) ,(4) 



where Ta is the nuclear thickness function, Ra is the 
nuclear radius and A the atomic number. 

We use the modified fragmentation functions 
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i/{Pn — Ai^n) is the rescaled momentum 
fraction, of hadron Hm originated from the fragmenting 



parton P„, z^ 



{L/X){hm/^En) is the rescaled mo- 



mentum fraction of the radiated gluon, AEn is the aver- 
age radiative parton energy loss and (L/X) is the average 
number of scatterings. The energy loss AEn is related to 



the gluon density of the produced medium via 
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and b is the impact parameter of the colhsion, r is the 
transverse plane location of the hard scattering where the 
partons are produced and h is the direction in which the 
produced hard parton travels in the medium, r and ip 
are the radial position and the angle that n makes with 
the radial direction, respectively. The function i(r, ip) 
can be interpreted as the medium influence time, that is 
the time during which the medium makes a parton loose 
energy, when it was produced at r and starts traveling 
through the medium making an angle ip with the radial 
direction. The parameter tq is the formation time for 
the medium gluons and po is the initial central gluon 
density, r is the time elapsed from the formation time 
and parametrizes the path length over the trajectory of 
the parton within the plasma. Since the gluon density 
profile is a rapidly falling function, the upper limit of 
integration can be safely set to infinity. Figure [T] shows 
t{r, ip) for a central collision. Notice that because of the 
dilution of the medium due to longitudinal expansion of 
the plasma, the biggest effect is on partons which have to 
travel less than the full length of the plasma, i.e. partons 



which are produced around two thirds from the center of 
plasma on the transverse plane. 

The average number of scatterings for a given parton 
is 
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where Aq is the parton mean free path. Since we want 
to consider the most central collisions, hereafter we set 
b — 0. The one dimensional energy loss {dEn/dL)id is 
parameterized as 
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The one-dimensional energy loss per unit length parame- 
ter, eg, is related to the mean free path Aq by epAo = 0.5 
GeV. We work with a value /zq = 1.5 GeV and refer the 
reader to Refs. [§-[l2[ for explicit details of the meaning 
and values of the introduced parameters. 

Figure [2] shows the average number of scatterings per 
unit length for a central collision. Shown is the case for 
£0 = 1-5 GeV/fm for which the maximum average num- 
ber of collisions is about 2.5. The average number of col- 
lisions grows with eo. As is clear from the figure, particles 
that experience the largest average number of collisions 
are the ones that arc produced around one third from the 
center and travel opposite to the radial direction, toward 
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FIG. 1: Medium influence time i(r, ip) for a central collision 
as a function of all the initial points and trajectories on the 
transverse plane for a central collision. The initial points are 
characterized by their radial position r and the angle that the 
parton trajectory makes with the outward radial direction ip. 
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FIG. 2: Average number of scatterings as a function of all the 
initial points and trajectories on the transverse plane for a 
central collision. The initial points are characterized by their 
radial position r and the angle that the parton trajectory 
makes with the outward radial direction (^. Shown is the case 
for eo = 1.5 GeV/fm. 



the interior of the medium. This effect can also be un- 
derstood as arising from the dilution of the medium due 
to longitudinal expansion of the plasma. 
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Given the above, the differential cross section for pro- 
duction of three hadrons with momenta hi, h2, h^ at 
mid-rapidity in A + A collisions is given by 
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where the new integrations with respect to the p -I- p case 
are performed over the overlap area of the nuclear colli- 
sion and over the angle ip that the direction of emission 
of the leading hadron makes with respect to the radial 
direction. We have correspondingly divided by the nu- 
clear overlap area ttR\ and by 27r in order to have the 
per nuclcon yield. 

The matrix elements Aiij^p^p^p^ in Eqs. ([1]) and ([TU]). 
representing the 2 — )- 3 LO QCD hard scattering am- 
plitudes at the parton level are divergent due both to 
coUincar and soft singularities. These kinematic di- 
vergences are universal and can be removed using well 
known techniques [l3|, |lj] . Here instead we apply angu- 



lar cuts to avoid the divergent regions. An alternative 
is to use an automated version of the dipole subtraction 
method [l^ in the MadGraph environment [1^1 ■ This 
avenue is being explored and will be reported elsewhere. 



Hereafter, A4i 



stands for the finite matrix ele- 



ments, after applying angular cuts. Throughout we use 
the factorization and renormalization scales /xj and fXr as 
Mr = M/ = 2 GeV. 
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FIG. 3: Cross section (da/dSf dSf ) for a symmetric hadron 
momentum configuration where 2 GeV < hi, /12 , /is < 3 GeV 
with y^sjvjv ~ 200 GeV. From left to right and top to bottom 
are the p -f p case and the A -f A cases for eo = 0.1, 1.5, 2 
GeV/fm, respectively. 



FIG. 4: Cross section {da/d62d.6f) for a symmetric hadron 
momentum configuration where 2 GeV < hi, /12, /is < 3 GeV 
with ^ysNN = 2.76 TeV. From left to right and top to bottom 
are the p 4- p case and the A -I- A cases for eo = 2, 3, 4 
GeV/fm, respectively. 
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FIG. 5: Cross section (da jdO^dO^) for an asymmetric hadron 
momentum configuration where 3 GeV < /ii < 4 GeV and 2 
GeV < /i2, /i3 < 3 GeV with ^/s]^ = 200 GeV. From left to 
right and top to bottom are the p + p case and the A + A 
cases for eo = 0.1, 1.5, 2 GeV/fm, respectively. 



FIG. 6: Gross section {da/d9^d6f) for an asymmetric hadron 
momentum configuration where 5 GeV < hi < 6 GeV and 2 
GeV <h2, /i3 < 3 GeV with ^/s]^ = 200 GeV. From left to 
right and top to bottom are the p + p case and the A + A 
cases for eo — 0.1, 1.5, 2 GeV/fm, respectively. 



III. THREE-HADRON CORRELATIONS 

Figure [3] shows the contribution of 2 — ?> 3 processes to 
the three-hadron correlation function. Shown is the cross 
section (dcr/dflf dflf ) for a symmetric hadron momentum 
configuration where 2 GeV < hi, /12, /13 < 3 GeV with 
a collision energy of -y/ijviv — 200 GeV, appropriate for 
RHIC energies. From left to right and top to bottom, 
the figure shows the p + p case and the A + A cases for 
values of the one-dimensional energy loss per unit length 
parameter eo = 0.1, 1.5, 2 GeV/fm, respectively. These 
values are chosen to correspond to the ones explored in 
Ref. H]. In these plots n/5 < 6^(6^) < 47r/5, 6tt/5 < 
^3(^2) — 97r/5. The excluded region corresponds, to 
events with only two hadrons in the final state, within 
our resolution. We also notice that when eo — > 0, the 
correlation function in the p + p case is recovered and 
that with increasing values of eg the intensity of the signal 
decreases, as expected. 

Figure m shows also the contribution of 2 —J- 3 processes 
to the three-hadron correlation function for a symmetric 
hadron momentum configuration 2 GeV < /ii , /12 , /13 < 3 
GeV but this time for a collision energy of y^SNN = 2.76 
TeV, appropriate for LHC energies. From left to right 
and top to bottom, the figure shows the p -I- p case and 
the A -I- A cases for values of the one-dimensional energy 
loss per unit length parameter cq = 2, 3, 4 GeV/fm, 
respectively. As in Fig.[3]for these plots 7r/5 < 62(0^) < 
47r/5, 67r/5 < 0^(9^) < 97r/5. We also notice from Fig.g] 



that with increasing values of eo the intensity of the signal 
decreases, as expected. 

Figures [5] and [6] show asymmetric hadron momenta 
configurations. Figure [5] shows the case where 3 GeV 

< hi <4 GeV and 2 GeV < /12, /13 < 3 GeV, whereas 
Fig. El is the case where 5 GeV < hi < 6 GcV and 2 GeV 

< /i2, h^ < 3 GeV, both calculated with ^/snn = 200 
GeV. Notice that the position of the two peaks on the 
away side remain at roughly 27r/3 and 47r/3 rad, as was 
the case for the symmetric configurations. Their inten- 
sity decreases with respect to the symmetric, lower mo- 
mentum case of Fig. |3l as the difference between away 
and leading particle momenta increases but the higher 
value becomes sharper. Also, the intensity decreases as 
the energy loss parameter increases. 

In order to extract information from the three-hadron 
correlation function in a two dimensional analysis, one 
possibility is to look at this object as a function of the 
angular difference Acp = 0f — ^2 of the away side particles 
for a range of angles of one of the away side-particles, say 
A9 = O2 ■ Figure [7] shows this correlation in a A -|- A en- 
vironment with y/SNN = 200 GeV for a leading hadron 
momentum 2.5 GeV < /ii < 4 GeV and away-side hadron 
momenta 1 GeV < /12, h^ < 2.5 GeV, integrated over a 
02 angular range 1.65 < A9 < 2.2 rad. Shown are the his- 
tograms obtained for eo = 1,2,3 GeV/fm normahzed to 
their cross section, ct, obtained by integration of the dif- 
ferential cross section over the above angular ranges and 
compared to preliminary data from PHENIX Q . Due to 
our angular cuts, meant to exclude collinear hadron pro- 
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FIG. 7: Cross section da/d9203 as a function of the angu- 
lar difference A(f) — 0^ — 02 for a leading hadron momentum 
2.5 GeV < /ii < 4 GeV and away-side hadron momenta 1 
GeV < /i2, hs < 2.5 GeV integrated over a A9 = ^^ angular 
range 1.65 < A6 < 2.2 rad. The histograms are normalized 
to their cross section, a, obtained by integration of the differ- 
ential cross section over the above angular ranges and corre- 
spond to three values of the energy loss parameter eo = 1, 2, 3 
GeV/fm. The calculation is compared to preliminary data 
from PHENIX. 



duction, the regions Ac/) ~ 0, tt rad are not accessible. 
However, as is clear from the figure, the angular region 
1.5 < A(f> < 2.7 rad is well described, particularly for the 
value of eo = 2 GcV/fm. 



cation of the peaks stays the same and their intensity 
decreases as the difference between away and leading par- 
ticle momenta increases. The intensity also drops as we 
increase the energy loss parameter. 

There are a few interesting questions that warrant fur- 
ther work, for instance, one could consider the contribu- 
tion of these 2 — ?> 3 processes to di-hadron correlations. 
This could happen when either one of the partons travels 
a large distance in the plasma and loses a large portion 
of its energy, so that its transverse momentum is out- 
side the transverse momentum window considered for the 
associate hadrons. Another possibility comes from the 
projection of the three-hadron angular correlation onto 
two-dimensions, i.e., when one integrates over one of the 
away side hadrons' angle. 

It would be interesting to compare these contribu- 
tions with those of genuine 2 — > 2 processes. Naively, 
one might expect that for the case when there are three 
hadrons in the final state, the partonic hard scattering 
must happen close to the center of the plasma in the 
transverse plane. On the other hand, one expects that 
contributions of 2 — >■ 3 processes to di-hadron correla- 
tions come from events where the partonic hard scatter- 
ing happens closer to the edge, so that at least one of 
the partons travels a large distance through the medium 
and has a high probability of being completely absorbed. 
A study of the transverse plane location dependence of 
these events may thus provide more information for the 
tomography of the Quark-Gluon Plasma than is possible 
with single or even genuine double inclusive hadron pro- 
duction. This work is in progress and will be reported 
elsewhere. 
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IV. CONCLUSION 

In this work we have calculated the three-hadron pro- 
duction cross sections in both p -I- p and A -I- A colhsions 
and investigated the effect of parton energy loss on the 
angular correlations between the produced hadrons. We 
have considered collision energies appropriate for RHIC 
and LHC. For the cases where the hadrons have same 
or not too different momenta, studied in this work, we 
observe two peaks on the away side, with respect to the 
leading hadron, at roughly 27r/3 and 47r/3 rad. The lo- 
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